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We study spin polarization in a split-gate quantum wire focussing on the effect of a realistic smooth 
potential due to remote donors. Electron interaction and spin effects are included within the density 
functional theory in the local spin density approximation. We find that depending on the electron 
density, the spin polarization exhibits qualitatively different features. For the case of relatively 
high electron density, when the Fermi energy Ef exceeds a characteristic strength of a long-range 
impurity potential Vdonors, the density spin polarization inside the wire is practically negligible and 
the wire conductance is spin-degenerate. When the density is decreased such that Ef approaches 
Vdonors, the electron density and conductance quickly become spin polarized. With further decrease 
of the density the electrons are trapped inside the lakes (droplets) formed by the impurity potential 
and the wire conductance approaches the pinch-off regime. We discuss the limitations of DFT-LSDA 
in this regime and compare the obtained results with available experimental data. 



Introduction 

The possibility of a spontaneous spin-polarization at 
low electron densities in low dimensional electron systems 
has attracted an enormous interest over the past years. 
The phenomena has been suggested to occur in various 
systems including quantum point contactsi (QPCs), two- 
dimensional electron gas {2DEG)^, quantum wire^ and 
open quantum dots^. Theoretical modeling with Hartree- 
Fock^, the density functional theory^i^i^ and Monte 
Carlo simulationsiS has reproduced low density spin po- 
larization in a number of systems. The mechanism driv- 
ing the polarization is the exchange energy dominating 
over the kinetic energy at low densities, making the spin- 
polarized state the energetically most favorable. The 
electron density necessary for this polarization to occur is 
generally very low, below ~ 3 x lO^'^m"^ (correspond- 
ing to the interaction parameter— = l/a^y^TmJ sa 3.2, 
where is the effective Bohr radius) as indicated in, 
e.g.fi^iH. At such low electron densities the electrostatic 
potential due to impurities in the donor layer can sig- 
nificantly affect the electronic and transport properties 
of the 2DEG. For example, Nixon et alJ^ showed that a 
monomode quantum wire is difficult to achieve because of 
the pinch-off due to a random potential from unscreened 
donors. This pinch-off is characterized by the critical 
electron density Uc , - the density where the 2DEG un- 
dergoes a metallic-insulator transition (MIT)iiii^ii^iii. 
The MIT causes a localization of the electron gas with 
an accompanying abrupt change in the conductanc o^^'^^ . 
Recent measurements of the thermodynamic magnetiza- 
tion in silicon 2DEG:iM^ found an enhancement of the 
spin susceptibility close to Uc- Interestingly, theoreti- 
cal consideration a^'^" indicate that the spin-susceptibility 
in a disordered potential increases for electron densities 
close to ric- These studies considered the general behavior 
of a 2DEG in an impurity potential but did not elaborate 
on specific geometries, e.g., quantum dots or wires. 

In the present paper we study spin polarization in a 
split-gate quantum wire focussing on the effect of a re- 



alistic smooth potential due to remote donors. For this 
purpose we, starting from a heterostructure and a gate 
layout, model GaAs/AlGaAs quantum wires within the 
density functional theory in the local spin density ap- 
proximation (DFT-LSDA) accounting for a long-range 
impurity potential due to ionized dopants. A gate volt- 
age applied to a top gate allow us to tune the electron 
density in the wire close to Uc. We find that depend- 
ing on the electron density, the spin polarization exhibits 
qualitatively different features. For the case of relatively 
high electron density, when the Fermi energy Ep exceeds 
a characteristic strength of a long-range impurity poten- 
tial Vdonors , the density spin polarization inside the wire 
is practically negligible and the wire conductance is spin- 
degenerate. When the density is decreased such that Ep 
approaches Vdonors, the electron density and conductance 
quickly become spin polarized. With further decrease 
of the density the electrons are trapped inside the lakes 
(droplets) formed by the impurity potential and the wire 
conductance approaches the pinch-off regime. 



Model 

We study the conductance and electron density of a 
GaAs/AlGaAs quantum wire with a realistic long-range 
impurity potential due to remote donors. In the wire, 
sketched in Fig. 1, the confinement is induced by two 
metallic side gates situated 700 nm apart on the top of 
the heterostructure. The heterostructure consists of the 
cap layer, the donor layer and the spacer. Electrons from 
the fully ionized donor layer form a two-dimensional elec- 
tron gas at a GaAs/AlGaAs interface situated at the dis- 
tance d2DEG below the surface. The confining potential 
from the donor layer is different in the leads and in the 
central region of the device. The leads, which extend to 
electron reservoirs at infinity, are considered ideal and the 
ionized donors in the lead regions are treated as a uniform 
layer with the density pd at the distance dd from the sur- 
face. Thus, the leads guides charges from the reservoirs 
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FIG. 1: (color online) Schematic view of the system studied. 
The heterostructure consists of (from bottom to top), a GaAs 
substrate, a 30-60nm AlGaAs-spacer, a 26nm donor layer and 
a 14nm cap layer. The side gates define a quantum wire and 
the top gate controls the electron density in the section with 
randomly distributed dopants Pd(r). 



to and from the middle region without any scattering. In 
the middle section the average donor density is still pd, 
but a gradual transition from the uniform donor density 
in the leads to random placement of the dopants is im- 
plemented. As a results, the electrons are scattered in 
the middle region of the wire due to the long-range im- 
purity potential. The donor potential on the depth of the 
2DEG is calculated according to the procedure outlined 
by Davies et ali^. An additional top gate in the middle 
section allows us to control the electron density in this re- 
gion. The potential on the top gate, Vg, ranges from zero 
voltage up to —0.07 V, which is the pinch off voltage for 
at least one of the spin species. The self-consistent po- 
tential at the GaAs/AlGaAs interface for different gate 
voltages is illustrated in Fig. [21 The shape of the im- 
purity potential depends on the spacer thickness. We 
consider two cases when the width of the spacer layer is 
respectively 30 nm and 60 nm. In both cases the aver- 

1 X lO^^m-^ 

rg — 0). This is achieved by choosing slightly dif- 
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1 X 1024m-3 and = 1.07 x lO^-^m-^ for re- 



spectively 30nm and 60nm spacers). The effective width 
of the wire in the 2DEG is 500 nm for both cases. For 
the case of 60nm spacer, the donors in the central region 
are situated further away from the 2DEG which results 
in a smoother profile in comparison to the SOnm case, c.f. 
Fig. [21(a), (c),(d) and (b),(d),(f). 

Using the Kohn-Sham formalism we write the Hamil- 
tonian for the quantum wire as21 



2m 



(1) 



where m* = 0.067me is the effective mass in GaAs, 
r = (cc, y) and a stands for spin up/down electrons (t, i). 
The total potential y (r) can be written as the sum of 
the classical Hartree potential, (r), the correlation and 
exchange potential, V^'^(r), and the external potential 
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FIG. 2: (color online) The self consistent spin-up/down po- 
tential along a slice in the middle of the wire (y = 0) for 
increasing gate voltages, Vg . (a), (c) and (e) show the 30nm- 
spacer sample while (b), (d) and (f) the 60nm-spacer. 



due gates, donors and Schottky barrier, Vext{^)- 

V^{r)^VH{r) + V,^,{v) + Ve.t{r) (2) 

With mirror charges at distance eg above the surface 
the Hartree potential is written as 



dr'n(r') | 
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(3) 



2DEG , 



where n(r) is the total (rt^ + n'^) electron density. Within 
the LSDA approximation the exchange and correlation 
potential is given by 
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For Exc the parametrization by Tanatar and Ceperly^^ 
was implemented. Finally, for Ve2,t(r) = Vgotes(r) + 
Vdonors{ic) + Vschottky{r) wc usc analytical expressions for 
Vgates (r)^ and Vdonors (i") (Rcfs^ audii respectively for 
the lead- and middle sections of the wire); the Schot- 
tky barrier Vschottky[Y) is set to 0.8cV. Using the recur- 
sive Green's function technique with mixed basis seii^^ 
we compute the conductance through the scattering re- 
gion (middle section) and the self-consistent electron den- 
sity in the system. Details of our implementation can 
be found in2i i26,27 g^j^j ^j-^^ procedure will only be briefly 
sketched here. The Hamiltonian Eq. ([T]) is discrctized on 



an equidistant grid and the retarded Green's function is 
defined as 
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(5) 



The electron density is integrated from the Green's func- 
tion (in the real space), 
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fpD being the Fermi-Dirac distribution. First we com- 
pute the self-consistent solution of equations (H])-® for 
an infinite homogenous wire by the technique described 
in2£. The converged solution for the infinite wire is used 
to find an approximation for the surface Green's function 
of the left and right leads. This approximation can be jus- 
tified because of a sufficient separation between the leads 
and the inhomogeneous potential in the middle section 
such that any inhomogeneous contribution to the poten- 
tial from the middle section is negligible at the leads. 
Next we apply the Dyson equation to couple the left and 
right surface Green's function and recursively compute 
the full Green's function for the middle section. We then 
iterate Eq. ([I])-® to find a self-consistent solution for the 
middle section. On each iteration step i the electron den- 
sity is updated from the input and output densities of the 
previous step, n™j(r) = (1 — e)n™(r) -|-en°"*j(r), e being 
a small number, ^ 0.05. Convergence is defined as a ra- 
tio between the relative change in input/output density 



at the iteration step i. 



n™) < 10- 



Finally the conductance is computed from the Landauer 
formula, which in the zero bias limit is 



G" = - — 



dET^iE) 



dE 



(7) 



where T'^ (E) is the transmission coefficient for spin chan- 
nel a. T'^{E) can be found from the Green's function 
between the leads^^. Calculations are done at zero mag- 
netic field. In order to find spin separated solutions a 
small magnetic field, ~ 0.05T, is applied for the first 
-100 (out of 1000-20000) iterations. Although the direct 
effect of the magnetic field is very small it is sufficient to 
lift spin degeneracy and for the converged solution to be 
spin polarized. The temperature in all simulations was 
chosen to IK. 



Results and Discussion 



The top panels in Fig. [3] show the spin up/down elec- 
tron densities and the spin polarization — at 
= 0. For this gate voltage the impurities in the middle 
section cause a clear modulation of the electron density 
with a negligible spin polarization inside the wire. How- 
ever, the electron density exhibits a pronounced spin po- 
larization near the wire edges. Because incoming states 
in the leads are spin degenerate, the polarization along 
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FIG. 3: (color online) Electron densities n\ (left and mid- 
dle columns) and spin polarization |n-^— (right column) for 
the 30nm-spacer sample at different gate voltages, Vg = OV, 
-0.03V, -0.07V (rows (a), (b), (c) respectively). 



the edges was quite unexpected. To understand the ori- 
gin of this spin polarization, we study an infinite ideal 
homogeneous wire, where the confinement is modeled by 
a parabolic potential. 



Vpariy) = Vb + y {uiyf , 



(8) 



Vq being the bottom of the parabola. Note that 
the parabolic confinement represents an excellent ap- 
proximation to the electrostatic potential from a gated 
structuresiii^iS^. At the same time, by changing the 
saddle point of the parabola, Vq, and the confinement 
strength, hui, it is convenient to control both the elec- 
tron density and the smoothness/steepness of the poten- 
tial. Self-consistent solution of Eqs. (U)-® can be spin 
degenerate and spin polarized (left and right panels in 
Fig. 0] (b), (c)). As for the case of quantum wire of Fig. 
[1] a small magnetic field was temporarily introduced to 
trigger spin degenerate solutions for some initial itera- 
tions. Figure m^c) shows a representative spin-polarized 
electron density in an ideal infinite quantum wire. As 
the electron density decreases at the edge of the wire, it 
becomes spin polarized and exhibits a spatial spin po- 
larization yielding a separation, d^ep, between the spin 
up/down densities. This is summarized in Fig. 3] for a 
series of wire configurations. Along the Vo-axis in Fig. 
mja) the width of the wire is held constant whereas the 
electron density grows as |Vb| is increased (we keep the 
Fermi energy Ep = 0). Conversely, changing the width of 
the wire along the if-axis by decreasing the confinement 
strength, huj, and keeping Vq constant, a more shallow 
wire is studied. The behavior of the spatial spin polar- 
ization presented in Fig. |3Ja) shows that d^ep increases 
as the electron density is decreased and the confinement 
becomes smoother. Note that this dependence of the dgep 
as a function of the electron density and the confinement 
strength is consistent with the corresponding behavior of 
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FIG. 4: (color online) (a) Spatial spin separation, dsep, at the 



degeneracy holds for energy levels far away from Ep, see 
Fig.Hb), (d). 

Having established that an infinite quantum wire can 
have a spin-polarized solution, we conclude that this so- 
lution is triggered in the finite quantum wire as well, even 
though the electrons injected into the middle part of the 
wire arc spin-degenerate (we stress that we always se- 
lect the spin-unpolarized solution in the leads). A word 
of caution is however in order concerning a reliability of 
the above predictions for the spatial spin separation near 
the wire edges obtained within the DFT-LSDA. Our re- 
cent comparison of the DFT-LSDA and the Hartree-Fock 
(HF) approaches demonstrates that the two methods pro- 
vide qualitatively (and in most cases quantitatively) sim- 
ilar results for electronic properties of ideal infinite quan- 
tum wires in the integer quantum Hall regime^. How- 
ever, in contrast to the HF approach, the DFT calcula- 
tions predict much larger spatial spin separation near the 
wire edge for low magnetic fields (when the compressible 
strips for spinless electrons are not formed yet). Note 
that a comparative study of two methods can not dis- 
tinguish which approach gives a correct result for d^ep 
for zero field. This question can be resolved by a com- 
parison to the exact results obtained by e.g. quantum 
Monte Carlo methods. We thus can not exclude that the 
predicted spin polarization near the wire boundaries as 
B = can be an artifact of the DFT-LSDA, and we defer 
this question to further studies. 

Let us now focus on the spin polarization in the cen- 
tral part of the wire. Depending on the electron density, 
we can identify three regimes with qualitatively different 
behavior. In the first regime the spin polarization of the 



boundary of the quantum wire vs the saddle point potential Vo electron density P„ 



n'\/(n' +n') and the spin po- 



and the wire width w. dsep is loosely defined as the distance 
between spin species at the level Q.hxn(y = 0). Vo and w 
can be ascribed as the electron density ant potential profile 
smoothness, respectively. In (b)-(f) the left and right panels 
show respectively the spin degenerate solutions (black lines) 
and spin-resolved solutions (red and blue lines), (b) and (d), 
electron density and band structure in the wire indicated (B) 
in panel (a), (c) and (f), electron density and band structure 
in the wire indicated © in panel (a). 



dsep near the edges of a quantum wire in perpendicular 
magnetic field^, where dsep also increases as the elec- 
tron density is decreased and the confinement becomes 
smoother. It should be noted that in the present case of 
zero magnetic field, d^ep shows a nonmonotonic depen- 
dence of the spin polarization and electron density/slope 
of confinement potential. This behavior of dsep is a man- 
ifestation of the subband structure in a quantum wire of 
a finite width. For example, if there is some energy level 
close to Ep, the exchange interaction effectively splits it 
into spin-up and spin-down subbands, see Fig. [lllc),(f). 
As a result, the total charge density profile shows a spa- 
tial spin separation dsep at the boundary of quantum 
wire (left part of (c) in figure 2]). In contrast, a spin- 



larization of the conductance Pq = |GT - GT|/(GT + G^) 
are negligible; for the 30nm spacer sample this is roughly 
between -0.025V< Vg < OV (row (a) in Fig. ^ 
while for the 60nm spacer sample this happens between 
— 0.045V< Vg < OV. In this regime the decreasing gate 
voltage causes a decreasing conductance (because of a 
reduction of the number of propagating subbands), but 
no significant spin polarization, except at the edges, oc- 
curs. The polarization at the edges is expected since 
the wire under consideration is wide (~500nm) and shal- 
low (minimum potential ~ — 5meV), which corresponds 
to the high polarization region of Fig. HJa). The self- 
consistent potential, shown for a slice along the middle 
of the wire in Fig. [U^a), (b) is well below the Fermi en- 
ergy, such that the characteristic potential fluctuations of 
the long-range impurity potential are much smaller than 
the average distance from the potential bottom to Ep. 

As the gate voltage becomes more negative, the wire 
undergoes a spin polarization in the central part, see Fig. 
I^b), (c). For the 30nm spacer sample (row (b) in figure 
[3]) this occurs for the gate voltages Vg < —0.025V and 
for the 60nm sample, Vg < —0.045V. The splitting re- 
sults in the fragmentation of spin- up/down densities into 
spin-polarized islands in the wire seen in row (b) of fig- 
ure [51 The onset of spin polarization is displaced towards 
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FIG. 5: (color online) (a) The average electron density, + 
n^, directly beneath the top gate. The arrow indicate the 
critical electron density, nc, for the MIT according to Eq. (|9]) 
for the 30nm spacer sample. For the 60nm spacer sample 
ric = 3.0 X lO^^m"^ is achieved for a gate voltage lower than 
-0.07 V. (b) The density spin polarization, P„; (c) the spin 
polarization of the conductance, Pg\ (d) the spin-resolved 
conductance C . 



a lower gate voltage for the thicker spacer. Fig. [IJb/c), 
but it is not clear whether there is any other qualita- 
tive difference between the two samples. To settle this 
would require further computations over more samples 
with varying spacer thickness and donor sheet configura- 
tions. However, because of the extensive computational 
efforts needed to find a convergent solution (sometimes 
requiring up to 20000 iterations) we were in position to 
study only two representative wires with spacers 30 nm 
and 60 nm. 

Finally we identify a third, non-conducting regime, 
corresponding to a metal-insulator transition (MIT); 
Vg < -0.04 V for the 30nm spacer and Vg < -0.07 V 
for the 60nm spacer. Conductance is pinched off and 
electrons are trapped in isolated pockets along the wire. 
Fig. [21(c)- This electron-droplet state has been analyzed 
thoroughly iiilii^i^^iM. Using arguments based on the 
screening of the impurity potential by the 2DEG, Efros 
et al^ gave an expression for the critical density, ric, 
where the metal-insulator transition (MIT) occurs. 
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(9) 



where is a numerical constant =0.11 and pd the donor 
density. We use this expression to find the approximate 
point for the MIT in our samples. For the SOnm spacer 
Eq. (O yields 5.9 x 

for the 60nm spacer Uc = 



IQi'^m-^ (see Fig. [5] (a)) and 



3.0 X lO^^'m-^ The predic- 
tions for Uc agree rather well with the numerical results 



that show that the wire undergoes a transition to a spin- 
polarized regime before the MIT occurs (see Fig. [5]) (b), 
(c). Thus, close to the pinch-off regime, the DFT-LSDA 
predicts formation of the spin-polarized electron lakes 
trapped in the minima of the long-range impurity po- 
tential. Such localized states might be relevant to the 
experimental observations of Bird at al. that provide 
evidence of bound state-mediated resonance interaction 
between the coupled quantum point contacts close to the 
pinch-off regime^. While microscopic origin on the effect 
is still under debate, our findings indicate that because 
the spin-polarized localized states trapped in minima of 
the impurity potential are generic feature of modulation- 
doped split-gate heterostructures, they might be relevant 
for the interpretation of the effect reported by Bird et 
am. 

Finally, a comment is in order concerning applicabil- 
ity of the method used. In our calculations we assume 
a constant chemical potential throughout the system. 
This condition is certainly violated in the regime close 
to the pinch off when the electrons arc trapped in iso- 
lated lakes containing an integer number of charges (i.e. 
in the Coulomb blockade regime of electron transport). 
It is now well recognized that the standard DFT-LSDA 
is not expected to work in the Coulomb blockade regime 
of weak coupling because of the spurious self-interaction 
errors caused by the lack of the derivative discontinuity 
of the exchange and correlation potentials in the stan- 
dard DFT— . The validity of the present method is lim- 
ited to the case of strong coupling when electron num- 
ber in the structure is not quantized (i.e. the Coulomb 
charging is unimportant) and the conductance of the sys- 
tems exceeds the conductance unit Go = 2e^//i (see Refi^ 
for a detailed discussion and further references). Be- 
cause of the uncorrected self-interaction of the standard 
DFT-LSDA; the electron lakes in the pinch-off regime do 
not contain an integer electron number and the calcu- 
lated conductance does not exhibit expected Coulomb- 
blockade peaks. We thus conclude that while the present 
calculations qualitatively capture the onset of the pinch- 
off regime, a quantitative description of this regime (ac- 
counting for the quantized electron number in the pock- 
ets as well reproducing the Coulomb blockade peaks in 
the conductance) would require methods that go beyond 
the standard DFT-LSDA scheme utilized in the present 
calculations. 

Let us now compare our findings with available ex- 
perimental results. Spontaneous spin-polarization at low 
electron densities has been probed in various systems, 
Refs^ii^iHiisas. Iii2^ Ghosh et al. studied the evolu- 
tion of the zero bias anomaly (ZBA)22 in 2DEGs for low 
and zero magnetic fields. The behavior of the ZBA was 
associated with different spin states in the 2DEG and 
measurements over different disorder configurations (cool 
downs) and temperatures indicated that the spin polar- 
ization observed is a generic effect for low density 2DEGs. 
The ZBA was most easily observed in a small disorder 
window between the metallic and insulating regime. This 
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is qualitatively consistent with the window of high spin 
polarization we find above (Fig. [5](b)-(d)). Further ex- 
periments on the ZBA in 2DEG»i^ suggested the forma- 
tion of localized magnetic moments due to spin polarized 
regions in the 2DEG as the electron density is lowered. 
This was understood as an effect of the potential due 
to background disorder which is similar to the impurity 
induced spin polarized droplets wc find in Fig. [3] (third 
column). Many of the observations ini^ were strongest 
for electron densities around 1 — 3 x lO^^m"^, a slightly 
lower electron density than wc find. In our case the av- 
erage density below the top gate for the onset of the spin 
polarization is 6.4xl0^''m~^ for the 30 nm spacer case 
and 5.1 X 10^''m^^ for the 60 nm spacer case. 

A direct measure of magnetization of the 2DEG at 
low electron densities were done i n^^i^^ for Si-Si02 het- 
erostructure 2DEGs. By modulating an in plane mag- 
netic field and measuring the minute current between 
gate and 2DEG the thermodynamic magnetization of the 
2DEG is found through Maxwell's equation a^^'^^ . Both 
Prus et ali^ and Shaskin et ali^ find that the spin 
susceptibility is critically enhanced prior to the metal- 
insulator transition in the 2DEG. It is, however, not clear 
from the experiment whether a spin polarized phase ac- 
tually exists between the metal and insulating phase or if 
there is only an increased magnetization in the metallic 
phase. 

Resonant inelastic light scattering measurements on 
GaAs single wells showed direct evidence for spin po- 
larization at low densities^^. Calculations using time- 
dependent local spin-density approximation in the same 
paper predicted a stable polarized state below an elec- 
tron density of 3.4 x lO^^'m"^. This is once again slightly 
lower than we encounter. 



Conclusions 



Using the spin density functional theory we have stud- 
ied spin-polarization of a 2DEG in split-gate quantum 
wires formed in modulation-doped GaAs heterostruc- 
tures focusing on the effect of the long-range impurity 
potential originating from the remote donors. We find 
that depending on the electron density, the spin polar- 
ization exhibits qualitatively different features in three 
different regimes. For the case of relatively high electron 
density, when the Fermi energy Ep exceeds a character- 
istic strength of a long-range impurity potential Vdonors , 
the density spin polarization inside the wire is practically 
negligible and the wire conductance is spin-degenerate. 
We find however a strong spin polarization near the wire 
boundaries. When the density is decreased such that 
Ep approaches Vdonors , the electron density and conduc- 
tance quickly become spin polarized. With further de- 
crease of the density the electrons are trapped inside the 
lakes (droplets) formed by the impurity potential and 
the wire conductance approaches the pinch-off regime. 
Experimentally, spin polarization prior to localization of 
the 2DEG has been suggested in^ii^iH. The electron 
density where we find spin polarization in the wire is 
roughly equal to what has been determined experimen- 
tally in GaAs/AlGaA o^^'^'^ . Direct measurements of the 
magnetization of the 2DEG in Si-Si02 heterostructures 
suggests an increased spin susceptibilitjsi^iii^ close to the 
MIT but it is not clear in these experiments whether an 
spin polarized phase, as we find, exists. 
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